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Abstract. - We examine one-dimensional quantum scattering of a Dirac particle off relativistic 
potential barriers. With proper considerations of Dirac sea, existence of anomalous tunneling at 
zero incident-energy is revealed for a particular type of relativistic potential having same magni- 
tudes and opposite signs for scalar and vector components. It is also shown that this leads to an 
exotic short range limit of the potentials. 



Since its conception, Dirac equation has been plagued 
OJoy difficulties related to the Dirac sea, and it is only after 
■j-Jthc proper taming of this monster of negative continuous 



spectra, that useful results are extracted from it. This 



^fact is demonstrated again recently by two example, one 
^Tbn the resolution of Klein paradox [1,2], and the other on 
'~^he resolution of Plesset's no-bound state problem [3-5]. 

In hght of the resurgent interest in the Dirac equation in 
^condensed matter physics [6,7], as well as of traditional 
lO interest in nuclear physics [8], we believe it is timely to 
report, in this Letter, a new result concerning a subtle 
effect of negative spectra on simple one ciinicnsional po- 
. tential scattering of Dirac particle off relativistic potential 
barriers. 



Here, we show that proper consideration of Dirac sea has 



intriguing ramifications for low energy scattering matrices. 
' In particular, it is shown to cause anomalous tunneling at 
zero incident energy for "5* = —V potentials", a particular 
type of relativistic potential having the same magnitudes 
H and opposite signs for scalar and vector components of the 



> 



potential. 

We further show, by considering the short-range limit 
of this potential scattering, that it leads to exotic low-pass 
gaussian wave filter whose non-relativistic kinematics limit 
is "delta-prime" point interaction, that causes discontinu- 
ity in quantum wave function itself but not in its deriva- 
tive [9] , which is distinct from conventional delta- function 
point interaction. 

We start by considering Dirac equation in one dimension 
that takes the following two-component form 
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Fig. 1: Examples of Dirac spectra. The graph on the left de- 
picts a case of m > — s, while the one on the right is of m < —s. 



where m, e stand for the mass and relativistic energy of a 
Dirac particle, and 5* and V are scalar and (time compo- 
nent of) vector potentials. We only treat time- symmetric 
systems, so the spatial components of vector potential is 
absent. The prime signifies the spatial derivative 

We first consider one dimensional potential barrier of 
constant height located in positive x region, formally given 

by 



V{x)=ve{x), S{x) = sQ{x), 



(2) 



where 6(x) is a Heaviside step function which is one for 
a; > and zero for a; < 0. We define "mass excluded" 
energy w by 



w 



m, 



(3) 



which we assume to be positive. The spectra of a Dirac 
particle inside the potential barrier is composed of two 
disjoint continuous spectra separated by a gap. Some ex- 
amples of the spectra are shown in FIG. 1. The scattering 
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wave functions at x < and x > are given, respectively 
as 
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with the free momentum k = \ — vn? and the momen- 
tum "p = \/{e — v)"^ — (m + s)^ inside the potential bar- 
rier. The expression (4) is valid for the case of continuous 
spectra for barrier region, e > \m-\- s\-\- v. For the case of 
|m -|- s| > e > — \m + s\+v, we need to make replacement, 
p — in with n = \/—{e — v)^ + (m + s)^. The case of 
e < — |m-|-s|-|-w corresponds to a particle under the Dirac 
sea, for which a p — > —p is needed in (4), but this case 
is soon to be sown irrelevant. Reflection and transmission 
rate is given by the squared absolute values of coefficients 
R and T, respectively, when p is real. When k is real, on 
the other hand, T is the amplitude of wave function at 
classically forbidden region. 

Expressing the momenta in terms of the energy w, we 
have 



k — \/ w{w + 2m), 

p = ^{w- 



- u) + 1m + s — u) . 
The matching of wave functions (4) at a:: = gives 
1 - i? = T, g(l + i?) = T 

with 



(5) 



(6) 



9 



w{w^1m-\-s — v) 1 w{w+2m+s — v) 



(w+2m){w — s — v) i y {w+2m)(s+v — w) 
Here, Q is the Giachetti-Sorace factor given by 
Q = 1 — Q{v — s — 2m — w)&{s + v — w). 



Q. (7) 



(8) 



that represents the exclusion of wave function to the bar- 
rier region x > when the energy w hits the negative 
energy spectra of Dirac equation with potentials s and 
V. It is technically obtained from the proper connection 
condition f(0-) — 0, x(O-) — constant, which is obtain- 
able as the n ^ CXI limit of {x — 1)" potential, for which 
f{0-) — 0, x(0_) — constant is found to be the cor- 
rect condition [5] . Fuller picture of this peculiar boundary 
condition may require the treatment of the problem with 
proper field theoretical setup as in [2], where the exclusion 
factor Q could be understood as a result of many-body 
Pauli blocking. 

The solution of the problem (6) is elementary, and we 
have 
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For large enough w, that satisfies the condition w > 
\m + s\ — m + V, the spectra inside the potential region 
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Fig. 2: Examples of Dirac spectra under scalar and vector po- 
tentials, s and V. Here, s is chosen to be proportional to v. The 
blue region represents the continuous particle spectra, while 
the red region is the "Dirac sea" spectra, which is occupied for 
the vacuum state. White region is the gap. Mass is set to be 
m — 1. Among the lines separating spectra, solid line repre- 
sents w — V + s, on which we have g = oo and the other line 
w = V — s — 2m on which we have g = 0. At the crossing point 

of these two lines, g takes the value g — ' 



is continuous, and we have patial transmission and reflec- 
tion specifled by (9) with (7). We naturally have unitar- 
ity relation |i?p + |Tp = l. As we decrease w down to 
the threshold energy w = v — m + \m + s\, p approaches 
zero, and g becomes cither zero (if s < —m and there- 
fore w = V — s — 2m) or infinity (if s > —to and therefore 
w = v + s). They respectively correspond to perfect reflec- 
tion R = 1 with Dirichlet boundary (p{0-) = or i? = — 1 
with Neumann boundary x(O-) = 0- 

Below this threshold, \m + s\—m + v > w > — |to-|-s| — 
m + V (or |e — v| < |to -I- s|, if it occurs with positive w), 
we have exponential wave function with decay constant 
K. The full reflection |i?| = 1 with quantum penetration 
< T < 2 to the classically forbiden area is observed. 
Note that there is no problem in having \T\ > 1 in this 
case, since the unitarity is guaranteed by decaying wave 
function e~'^^. At the "Dirac sea" threshold, w = v — m — 
|to -|- s|, K approaches zero, and g becomes either inflnity 
(if s < —TO and therefore w = v + s) ot zero (if s > —m 
and therefore w = v — s — 2m). They again respectively 
correspond to perfect reflection R = —1 with Neumann 
boundary x(O-) = or i? = 1 with Dirichlet boundary 
V.(0_)=0. 

Below the Dirac sea threshold, w < —\m -\- w\ — m + v 
we have perfect reflection with Dirichlet condition R = \ 
as a result of Giachetti-Sorace factor. 

When the energy w approaches 0, we have k — > 0, that 
signifies the quantum penetration length to classically for- 
bidden region x > becoming infinite. However, (7) tells 
us that we have g = Q, and thus no penetration amplitude 
T — Q and, as a result, the perfect "classical" reflection 
i?= 1. 

Above statements are true in generic case, depicted in 
the left hand graph in FIG. 2, for example, but there is an 
exception to the case when potentials v and s are related 
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Fig. 3: Schematic representations of relativistic potential bar- 
rier scattering 

hy s + V = 0, depicted in the right hand graph of FIG. 2. 
In this case, there is a cancelation in the expression for g, 
(7) and we have 



w+ 2m— 2v 
w+ 2m 



2v— w— 2m 



2m 



(10) 



which is finite as a results of "merging" oi w = v + s line, 
on which g = oo holds, and w = line, on which g = 
holds. 

This means that for this special case of opposite-sign but 
equal-magnitude scalar and vector potential, s -|- = 0, 
we have singular infinite-range penetration limit k — > 
with finite amplitude < |T| < 2 for zero energy barrier 
reflection w — > 0. We rush to note that this poses no 
paradox of any sort, since wc still have full reflection \R\ = 
1 albeit with some nontrivial phase for R. This is a subtle 
but an exotic exception, nonetheless, whose significance 
soon becomes obvious in the following. 

We now consider one dimensional scattering by square 
well of constant height potential with spatial extension L 
(FIG. 3), formally given by 

V{x) = ve{x)e{L-x), S{x) = sQ{x)Q{L-x). (11) 

The scattering wave functions ata;<0, 0<a;<Z/ and 
L < X are given, respectively as 



1 

ik 



R 



= T 



m+e+s— i; 

1 

ik 



1 

— ik 

B 



-ikx 



—ip 



-ipx 



(12) 



In a similar manner to the previous case, smooth ccmnec- 
tion conditions for both large and small components at 
a; = and x = L gives 



1-R=A-B 
g{l + R)=A + B 

Ae^P^ - Be-'P^ = 



rji^ikL 



(13) 



Ae'P^ + Be-'P^ = gTe'''^ . 



Elementary calculation yields the following expressions 
for transmission and refiection amplitudes; 



il.'L 



COS pL — ^{1/g + g) sinpL' 



cos pL — + g) sin pL 



(14) 



This expression is literally valid for the energy w > |m -I- 
s\ + V. For the energy \m + s\ + v > w > —\m -I- s| -I- v, 
we have to make replacement p = in as before, which will 
result in the replacements cos pL cosh/tL and sin pL — » 
isinh kL in (14). This expression is reduced to T = 0, 
R = I for the energy w < —\m + s\ + v with which wc hit 
the Dirac sea spectra inside the potential barrier, where 
wc have Q = 0, thus .9 = 0. 

Wc look at the low energy limit of the scattering matrix 
T. Gcnerically, for the case of s ^ —v, we have the quan- 
tity g that approaches to zero as we take w — > limit, 
causing the divergence of 1/g, which guarantees the per- 
fect reflection 



T^O, R^ 1 as 



(s + v^O). 



(15) 



This simply is a exact expression of the intuitive statement 
that a generic obstacle works as a reflecting block for low 
energy projectile, or in other word, if we hit any barrier 
too slowly, we are bound to get reflected all the time. 

However, for the special case of s + v = 0, g takes the 
form (10) after cancelation of w in both denominator and 
enumerator of (7), and we have g = flnite and k (or 
p ^ 0) as we take w — > limit. We therefore obtain, from 
(14), a peculiar limit 



1, R^O as w^O {s + v = 0), 



(16) 



which signifies an anomalous full transmission at zero en- 
ergy. This is particularly intriguing for the case of de- 
caying wave in the gap region, in which k is real, where 
decaying length 1/k becomes infinity at w limit. 

The situation is immediately understood by inspecting 
the illustrations in FIG. 4, Here, the graph in the left 
depicts a generic case that has normal perfect refiection 
a,t w limit, while the graph in the right shows the 
anomalous zero energy transparency. The reason behind 
this transparency lies in the enhanced long range tunnel- 
ing inside the barrier, which occurs because, at w ^ 0, 
the energy approaches to the threshold of negative contin- 
uous spectra that exists right below w = for s + v = 
potentials. The presence of Dirac sea not only induces the 
perfect reflection for w > with v > \s\, for example, it 
also affects the decaying length and induces the anoma- 
lous tunneling, and transmission at w ^ limit for the 
case of s -|- u = potentials. 

Let us now consider the L — > limit of relativistic scat- 
tering. Straightforward limit will, of course, lead to dis- 
appearance of the barrier, T — > 1. Limit L ^ with 
constant volume integrals. 



V s 

on the other hand, leads us to 
T = 



(i-0). 



(17) 



cosp+ ^smp[{v+s)/K+ {v-s)K] 
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the Dirac equation (1) by eliminating small component x 
in the form 



d 1 d 
dx 2m* dx 



tp + Ulp — voif, 



(24) 



with effective mass to* and potential U defined by 



Fig. 4: Transmission rate as function of Idnetic energy. The 
graph on the left is for normal case, and the one on the right 
is anomalous case, s = —v. 



^sm (i[{v+-s)/K~{v^-s)K] ,„ 

n. = ; ; : — TT^^ : r— rH(— W— |s|) 



C0S/3+ ^ sin/3[(u+ s)/K+ (v- s)K] 



+Q{v+\s\), 



(18) 



with B = Jv^ — s^— i\/ — and K = . , , „ . 

We can interpret this result in terms of relativistic point 
interactions specified by the boundary condition which is 
a most general time-reversal symmetric one [10] 



\x{0+) W a \x{0^) 



(19) 



with — u^u^ — 1. The scattering off the point interac- 
tion (f9) is given 



, (20) 



which allows the identifications 



_sSin/3 _ sin/3 

u+ = {s + v)^^, u_ = (s-w)— (21) 

The special cases s ~ v and s — —v can be considered as 
the limiting cases of (18), and we have, for s — v 



1 iv/K 
T=- — , R= ' 



l + iv/K' ' l + iv/K' 
while, for s = — u, we have 



„ 1 „ —ivK 

T^- — , i? = 



l + ivK' 



I + ivK 



(22) 



(23) 



If we take the non-relativistic limit in kinematics, K 
k/{2m), these two cases are exactly identical to the scat- 
tering form delta and delta-prime point interactions [9], 
which represent high-pass and low-pass wave filters, re- 
spectively. Note that constructing non-standard point 
interactions, that results in (23), within non-relativistic 
framework involves highly singular procedures [11,12]. 

Finally, we ask a question whether we can construct an 
analogue of the phenomena we have found in the frame- 
work of non-relativistic Schrodinger equation. We rewrite 



w 

TO+- 



u^s + v. 



(25) 



Assuming the conditions w << m and \S ~ V\ << to, we 
obtain Schrodinger equation with effective potential which 
is given by the sum of vector and scalar potentials. This 
is nothing but the true non-relativistic limit. However, 
we obtain non-standard low-energy limit by assuming the 
non-relativistic kinematics w << to in conjunction with 
strong relativistic potentials [S* — ^ to. Specifically, we 
can reproduce anomalous transmission from Schrodinger 
equation (24) by setting S — —V which results in to* « 
m ~ V and U — 0. This means that we can construct 
a purely non-relativistic model of anomalous scattering 
and delta-prime point interaction with just effective mass 
and no potential. Readers are warned, however, that this 
non-relativistic analogue scheme works only to an extent: 
When 5* — F is negative in sign and so large, we obtain 
negative value for the effective mass to*. For this bona fide 
relativistic dynamics, non-relativistic analogue (24) does 
not make sense, and therefore does not exist. 

We have shown that a set of relativistic potentials hav- 
ing the property of same strength but opposite signs for 
scalar and vector components displays anomalous full tun- 
neling and transparency at zero energy, while barrier starts 
functioning at higher energy. The short range limit of this 
phenomenon leads to a smooth relativistic realization of an 
exotic point interaction, delta-prime, that conventionally 
requires singular and esoteric constructions within non- 
relativistic dynamics. It has been pointed out [13] that 
in three dimensions, the "S = — relativistic potentials 
have an esoteric property called pseudospin symmetry [14] 
that has been found to play important role in the degen- 
eracy structure of nuclear levels. Current work shows that 
there is yet another aspect to this pseudospin symmetric 
limit of relativistic potentials, which is revealed only in 
one dimensional systems. 
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